PIONIC CONTRIBUTION TO NEUTRINOLESS DOUBLE BETA 

DECAY 



J. D. VergadosW<( 2 \ Amand Faessler^ 3 ), and H. TokiW 
^'Physics Department, University of Ioannina, Ioannina, GR 451 10, Greece, 
™ Theory Division, CERN, Geneva, Switzwerland 
' 3 - ) Institute fur Theoterische Physik,Universitat Tubingen, Germany, and 
' 4 ' RCNP, Osaka University, Osaka, 567-0047, Japan 
(Dated: November 11, 2009) 

It is well known that neutrinoless double decay is going to play a crucial role in settling the 
neutrino properties, which cannot be extracted from the neutrino oscillation data. It is, in particular, 
expected to settle the absolute scale of neutrino mass and determine whether the neutrinos are 
Majorana particles, i.e. they coincide with their own antiparticles. In order to extract the average 
neutrino mass from the data one must be able to estimate the contribution all possible high mass 
intermediate particles. The latter, which occur in practically all extensions of the standard model, 
can, in principle, be differentiated from the usual mass term, if data from various targets are 
available. One, however, must first be able reliably calculate the corresponding nuclear matrix 
elements. Such calculations are extremely difficult since the effective transition operators are very 
short ranged. For such operators processes like pionic contributions, which are usually negligible, 
turn out to be dominant. We study such an effect in a non relativistic quark model for the pion and 
the nucleon. 

PACS numbers: 12.60Jv, 11.30Er, 11.30Fs, 23.40Bw 



INTRODUCTION 



The discovery of neutrino oscillations can be considered as one of the greatest triumphs of modern 
physics. It began with atmospheric neutrino oscillations [l| interpreted as — ► v r oscillations, as 
well as v e disappearance in solar neutrinos 0. These results have been recently confirmed by the 
KamLAND experiment which exhibits evidence for reactor antineutrino disappearance. As a 
result of these experiments we have a pretty good idea of the neutrino mixing matrix and of the 
two independent quantities Am 2 , e.g m\ — m\ and m\ — m\. Fortunately these two Am 2 values are 
vastly different, 

\Am 2 21 \ = \m\-m\\ = (5.0 - 7.5) x l(T 5 (eF) 2 

and 

|Am 2 2 | = \m\ - m\\ = 2.5 x lCT 3 (eV~) 2 . 

This means that the relevant L/E parameters are very different. Thus for a given energy the 
experimental results can approximately be described as two generation oscillations. For an accurate 
description of the data, however, a three generation analysis ^ s necessary. 

We thus know that the neutrinos are massive, with two non zero Am 2 , and they are admixed. 
We do not know, however, whether they are Majorana, i.e. the mass eigenstates coincide with their 
antiparticles, or of Dirac type, i.e. the mass eigenstates do not coincide with their antiparicles. 
Furthermore we do not know the absolute mass scale as well as the sign of Am§ 2 ■ The first question 
can be settled by neutrinoless double beta decay (Ov (3/3— decay). The second will also, most likely, 
be settled by this process. 

We should stress, of course, the fact that the light neutrino mediated process is not the only 
mechanism available for Qv(3f3 [6(. Among those are some which involve heavy intermediate particles. 
These lead to very short ranged two body effective transition operators, which must be dealt with 
care, due to the presence of the nuclear hard core. To this end three treatments have been proposed: 

• Treat the nucleons as composite particles (two nucleon mode). 
This can be done in the context of non relativistic quark model or simply by assigning to the 
nucleon a suitable form factor [7| ■ 



Consider the possibility of six quark cluster in the nucleus 



• Consider other particles in the nuclear soup. 
The most prominent are pions in flight between the two interacting nucleons @ 

In the present study we will examine the last possibility. This was examined long time ago 0] and it 
was revived in the context of R-parity violating supersymmetry a decade later (lot HH , as well as 
recently [l3j|. It was shown that in the context of R-parity violating supersymmetry the pion mode 
is more important than the two nucleon mechanism. The same conclusion was reached recently in 
the context of effective field theory [14j. 

In the above treatments the pions were treated as elementary particles. This approach is reason- 
able in particle physics, but one knows, of course, that the hadrons involved are not elementary. 
Furthermore a crucial factorization approximation has to be made, by inserting only the vacuum as 
intermediate state, (see Eqs (|82|) and (|85|) below). Finally, even though the hadrons are elementary, 
in the interesting case of the pseudoscalar coupling an assumption had to be made about the quark 
mass, taken to be the current quark mass. 

In this work we are going to adopt a different procedure. The hadrons will be assumed to have a 
quark substructure in the context of the harmonic oscillator. In the harmonic oscillator approxi- 
mation the internal degrees of freedom can be separated from the center of mass motion. In this 
approach one derives the effective operator at the quark level by a suitable non relativistic expansion 
of the elementary amplitude. In some processes in our formalism one extra qq pair must be produced. 
This can can be achieved either through the weak interaction itself or via the strong interaction. 
The net result is that, in this new approach, one obtains new types of operators, including some 
that are non local at the nucleon level. One must weigh these advantages, however, against possible 
shortcomings of the need for a non relativistic reduction of the transition operator at the quark level. 



THE CONTRIBUTION OF PIONS IN FLIGHT BETWEEN NUCLEONS 

As we have mentioned in the introduction when the intermediate fermion, e.g. the Majorana 
neutrino, is very heavy the transition operator becomes very short ranged. In this case the usual 
two nucleon mechanism may be suppressed due to the nuclear hard core and the contribution of 
other particles in the nuclear soup, such as pions, may dominate. These mechanisms at the nucleon 
level are illustrated in Fig. [TJ 

The two body double beta decay operator, associated with heavy intermediate particle exchange, 
will be normalized in a way which is consistent with the light intermediate neutrino. We begin with 
the intermediate heavy neutrino. Then: 

r\ v — ^Vn S{r 1 -r 2 ) (1) 

r m e m p 

L, R stand for leftt handed and right habded currents respectively with 

-< m v y L R Trip 

Vv = , V N =^ — >- ( 2 ) 

m e rriN 

The corresponding expression in momentum space becomes: 

r/ u — 6{Ti - rj^ra - r 2 ) <^ r/ N ' llpp (3) 

r m e m p 



Mpp = -— ^^(Pl + P2 Pi - P2)^(Pl,P2,Pl,P2) (4) 

[Zn ) 

The function A(pi, P2, pi, p' 2 ) depends on the assumed mechanism for the neutrinoless double beta 
decay. 




FIG. 1: The double beta decay of two neutrons into two protons at the two nucleon level (a) arising when 
all the intermediate particles at the quark level are very heavy. The double beta decay of a neutron with 
the simultaneous production of a n + , which is then absorbed by another neutron converting it into a proton 
(b) (one pion mode). A neutron can also be converted into a proton and a tv~ . The ir~ then double beta 
decays into a ir + , which subsequently is absorbed by another neutron converting into a proton (c) (two pion 
mode). 



The factor V^ R is not usually included in the nuclear matrix element. The factor " R ^" p will be 
absorbed into the effective nuclear operator, while the factor ^ will eventually be included in the 
effective coupling, as will be discussed in this work. 

With the above expressions the formula for the life time due to heavy intermediate neutrinos in 
left handed V-A theories can be cast in the form: 



[if 



/2J 



a 



01 



(5) 



The two nucleon contribution ( jr^) 2 Q,p — Qqt was inserted in the above equation merely for com- 
parison. 

The case of other heavy intermediate particles, as those encountered in the R-parity violating 
supersymmctry can be handled in a similar fashion: 
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In both cases: 



=Y,T + (i)T + (j)a i . ( r j ^F[ k \x v ) (8) 



i<j 



r 



M!?=^T+(i)T + (j)[3a i .hj<Tjfij-<ri-<Tj] ^i^K) (9) 
i<3 

Where Rq is the nuclear radius, x^ = m^rij and 

F^\x) = e- x ,F 2 (1) (x) = (x 2 + 3x + 3)e~ x , f[ 2) (x) = (x - 2)e~ x , F 2 (2) (x) = (x + l)e~ x (10) 
The function A(p±, P2, p' 1; p' 2 ) depends on the pion mode under consideration. 

THE 2-PION MODE 

The spin dependence of the transition operator is in this case trivial. So we will focus on the 
orbital structure of of the operator The function A(p%, p2, p[, p' 2 ) is independent of the momenta 
in the standard V-A theory as well as in the case of the scalar (S-S) theory. It is, however, a 
model dependent function in the case of psedoscalar (P-P) interaction encountered, e.g., in R-Parity 
violating SUSY mediated double beta decay. In the last case we find 

A = -~A X .A 2 (11) 

where Aj is the amplitude resulting from the non relativistic reduction of the pseudoscalar involved 
in the d — > u coupling, i.e. 

U{j)'i)^d{pi) -> A,;. <Ti (12) 

where <Xj is the spin of the quark i and 

M = -^-Pi - ^-P'i (13) 
We find it convenient to rewrite them as follows: 



2 ~~ V2 2m d ~ 2mJ + 2m d ~ 2m u '~2 ' ' 

Where q = P T is the momentum of the pion in flight between the two nucleons and p and p are the 
relative internal momenta (see next subsection). One normally ignores at this level the momentum 
carried away by the two leptons. The 2ir 0v — [3(3 decay contribution in the case of heavy Majorana 
neutrino or any other Majorana fermion is explicitly shown in Fig. [2] 

Orbital integrals in the two pion exchange. 

The pion wave function is given by: 

1/2 

^(Q,p) = v^(2V2j (27r) 3/2 5(V2Q-P,)<Mp) (16) 




FIG. 2: The 0vf3/3 decay of pions in flight (2ir mode of Fig. [1} illustrated at the quark level. / stands for a 
effective exchange of a heavy Majorana fermion ( heavy neutrino or, as in R-parity violating supersymmetry, 
a neutralino, gluino etc). The ellipses merely indicate that the pion is a bound state of two quarks. 



where P w is the pion momentum and 

P = A= (P2 - Pi) , Q = \ (P2 + Pi) (17) 

with pi and P2 being the momenta of the quark and antiquark participating in the pion. This wave 
function is normalized in the usual way: 

-< ^pJV>p'„ >-= 2£ 7r (2 7 r) 3 <5(P 7r - p;) (18) 
4>n{ p) is described by an Is harmonic oscillator state. In momentum space it takes the form: 



^(p) = ^(0)e-^ 2 )/ 2 , <j> v (Q) = J (19) 



Thus the orbital matrix element in this case takes the form: 



ME 2n = M^)H{^ - P;) , M 2 „ = 1 2 %* fg{x), f£(x) = ^ ( 2 °) 

Z7TV Z7T Oj\r 



where x = b n and are the harmonic oscillator (HO) size parameters for the pion and the 
nucleon respectively. We have decided to introduce the ratio x as a variable to be adjusted. 
In V-A theories after incorporating the spin we find: 

-^M 27 r = c 27T ml (21) 



with 



^ = w,<^,^ 1 -" 1 "^ fS(x) (22) 



where -< |o"i.<t 2 | >-= —3 is the spin ME. One now can construct the effective transition operator 
in coordinate space at the nuclear level. The effective coupling in V-A theory is given Q by: 

2 1 1 1 



Or 



2m N J 47r 6m£ f£ 



2 

a 2 n = xtj fl N N c ^ (24) 

J A 

Using f 2 NN = 0.08 and 6jv = 1-0 fm we find a 2n = 0.013 and 0.11 for x = 1.0 and 0.5 respectively. 
For the scalar interaction one gets the value /J/4 with the value of f$ depending on the specific 
particle model. 

The dependence of the results on the pion size parameter is exhibited in Figs [3] 
In the case of the pseudoscalar coupling, since the pion has spin zero, we encounter the combination: 

(Ai. <r 1 )(A 2 . <r 2 ) =► -i(Ai.A 2 )( ax. <r 2 ) (25) 
In this case one can show that the orbital amplitude is 

= Q(«3 + «S) ~ \{** -ufbW) £\x), = ± (26) 

Where q the momentum of the propagating pion and 

Kd = o — r~ ' Ku = n — ~T~ ( 27 ) 

The above equation can be rewritten in a way that the pion propagator is manifest: 
1 2m T 



M 2jI = 



27rV27r bl 



N 



In other words there appear two terms and c 27r &^(<7 + m 2 ,) with 



C - = V^K^J^ W d + Kl) + ^ - ^ fhlml ) (29) 

4, = -^^4^5^ " ^ a «/£(*) (30) 

The first gives rise to an effective operator similar to that of the V-A theory with a coupling 

2 1 f2 1 tWf^f 1 /^,^^, 1 ^ n2i,2 2 



with -< |<Ti.<T2| >-= —3. 

The second term, contributing when the u and d quarks are not degenerate, yields a coupling 
a27r(^i7r) where: 

a 27T (n ln ) = -^i/^^i/P)^)!^ _ Km)2) x | CTi . CT2 | ^ (32) 

J A V 27T ™p JTlpON 6 

which is associated with the operator with one pion propagator less, i.e. that encountered in the l7r 
mode (see below). Such an operator is absent in the elementary particle treatment, even though the 
quarks are assumed to be non degenarate. 




THE 1-PION MODE 



In this case a positively charged pion, produced in virtual double beta decay of a neutron into 
a proton, is absorbed by another neutron converting it into a proton. At the quark level the first 
of these steps is exhibited in Figs 0JH1 In these figures a qq pair is created out of the vacuum. In 
the first two figures this is achieved as, e.g ., in a gluon exchange [TEj or a multigluon exchange 
simulated in the 3 Pq model 17|.[l8l].[l9]. The latter is a fairly old model, which still continues 
to be successfully applied in the description of meson decays [2(j. In Fig. [6] this pair is created by 
the weak interaction itself. 
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FIG. 4: The pion mediated Qu j3(3 decay in the so-called In mode. At the top we show the diagram in which 
the quarks of the pion are spectators , i.e. the heavy intermediate heavy fermion f is exchanged between the 
other two quarks, x indicates that a qq pair is created out of the vacuum in the context of a multigluon 
exchange. We will call it direct diagram. 



The orbital part at the quark level 

Orbital wave functions in momentum space are expressed in terms of Jacobi coordinates: 

1/2 

ifa* = (2V2) (2^) 3/2 S - P ff ) M p) (33) 



1/2 

^p=(3V3) (2^) 3/2 ,5(V3Q-P) (j>{ £, 77) (34) 
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FIG. 5: The same as in Fig. [4] involving the exchange diagram. In this case the quark involved in the pion 
participates in the exchange of the heavy fermion f , co-operating this way with another quark belonging in 
the nucleon. 
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FIG. 6: The same as in Fig. [5] but in a novel mechanism, i.e. one in which the qq pair is produced by the 
weak interaction itself. 



1/2 

%h' = (3V3) (2tt) 3/2 5 (V3Q' - P') <H 17') 



Where P^, P and P' are the momenta of the pion and the two nucleons respectively and 



(35) 



V2 



(Pi - P2) , V 



V6 



Pi + p 2 - 2p 3 ) , Q 



V3 



(Pi +P2 +P3) 



(36) 



r = 4(p't - P2) , »?' - 4=(p'i + P2 - 2 P D , Q' = -^(pi + P 2 + PD (37) 



P = ^(P3 - P4) > Q* = ~f= (P3 + P4) 



(38) 



Where pi, i = 1, 3 are the momenta of the three quarks of one nucleon, p' l5 p' 2 , p' 4 the momenta of 
the three quarks of the other nucleon and Pi,p' 3 are those of the quarks involved in the pion. This 
notation was chosen since the interaction preserves the fermion lines Pi < — ► p' 
The above wave functions were normalized in the usual way: 



-< 0p|v*» >= (2tt) 3 s(p - p') , -< ihJM y= 2^(27t) 3 <5(p t - p;) 



(39) 



The internal wave functions are given by: 



0(0 = 0(O)e"^« 2 )/ 2 , m = 4/-%= etc (40) 

The pion wave function has already been defined (see Eq. (fT9| ). except that sometimes we will 
write: 

<M0) = <K0):r 3 / 2 , x = ^ (41) 

OAT 

The integrals over the momentum variables Q, Q and can be trivially performed due to the 5 
functions. Thus the orbital integral becomes: 

I m = (2tt) 3 S(P - P' - P^M (42) 

M = ^7^^ J d 3 £d 3 £d 3 V d 3 rj'd 3 p4>( i, r,)4>( r/O-M p)n p0 (43) 

where Qpp depends on the mechanism involved as we now discuss. 

1. The qq pair is created by the 0^/3/3 operator ( Ou qq case) 

The case in which the qq pair is created by the 0vf3(3 operator (see fig. [H]). Then up to terms 
linear in the momentum the effective operator takes the form: 

ujsiv) = 04- I 7T— + , ^ ' ) (scalar and vector) (44) 



up = 0-3. ( ^- 3 — | (pseudoscalar ) (45) 

\ 2m u 2m d / 



uja = i((T3 x cr 4 ). ( — ^- ) (Axial current) (46) 

V 2m u 2771^ I 

It is, of course, understood that the scalar and pseudoscalar must be multiplied by suitable 
coupling constants. The full operator takes the form: 

ttpp-K = 7^3^ (l>a -Pa ~P4 -Pij 5 (pi ~Pi) $ (p 2 - p 2 ) j,-5,F,P,^4 (47) 
The product of the three S functions can be cast in the form 

6 ( P - P - R 



'tt) $ - p'-A 6 (p 2 - P2) = 

S(P-P'- P.) <5 (V2(€ - C')) ^ f-^fo -»/) + §) 



By setting £ = £ and 77 = 77 + \ / 4q we get 



/ 4 (3 (g- V2p) m d + m u (~5q - 2V677 + 2p N )) \ 

VS(V) ="4- (48) 

\ bmdm u / 



^ 6m d m u J 

( ai (3 (q - V%p) m d + m u (5q + 2VEr) - 2p N ))\ ..... . 

w A = z(cr 3 x er 4 ). ^ ^ i (Axial) (50) 

V 6m d m u J 

After the integration (see next section) we get: 

f (3qm d + m u (2p N - 3q))\ 04-PiV 

ws(y) = cr4- 7 =^ ( 51 ) 

V Qm d m u ) m p 

/(3qm u - m d (3q+2pjv))\ ofc-PiV 

= o- 3 . (52) 

V om d m u ) m p 

■ ( v \ f (3q(m d + m u ) -2m u p N ) \ -3q + pjy 

UA = H°3 X (7 4 ). =>■ l(cr 3 X (7a). (53) 

\ bm d m u J m p 

The last expressions result in the case of the constituent mass for the quarks, m u = m d = m p /3. 
In the above equations: 

P + P' 

PN = , q = P - P' = P w (54) 

2. Double beta decay and strong qq production ( 3 Po qq case). 

In this case one needs the collaborative effect of the Qvj3j3 interaction acting between quarks 
together the strong interaction, which creates a pion out of the vacuum (a' la 3Pq model or 
multigluon exchange) : 

= . 9 ; cr 4 .B <5(p 4 + p 4 ) , B = p 4 -p 4 (55) 

where g' r a dimensionless constant proportional to the parameter g r = 13.4 ± 0.1, which is 
known from experiment. One finds 

, _ 3^(2^ + 3) 3/2 
9r ~ 5r 80v / 2^(0)7r 3 / 2 m P V^ (56) 

where 5where 0^(0) is the pion wave function at the origin. 

• The direct term in the one pion contribution. 
In this case (see fig. |4j none of the two interacting quarks participates in the pion as 
defined above. Thus we get: 

^/3/37r = J^3 S (Pi +P2 -Pi ~ P 2 ) (P3 ~ P 4 ) S (p 4 + p 4 ) CT 4 . (p 4 + p 4 ) (57) 

The product of the above three S functions can be cast in the form 

S (P-P - Pn) S (pi+ P 2 -Pi -P 2 ) S (p 4 + P 4 ) 

The first of these 5-functions expresses momentum conservation. Going into the Jacobi 
variables we find: 



(2tt) 3 



.f2q-2V&V-^V2p + 2p N \ Aq + 2V&ij - 3^2p - 2p N 

5 I g U. g (58) 



We find it convenient to use the above d functions to obtain: 

2q + 3V2p - 2p N , 2q-3V2p + 2 PN 

v = Ul ' v = tVE (59) 

One finds: 

<x 4 .B= <t 4 . (q-V2p) (60) 

Furthermore A-terms, appearing in the case of the pseudoscalar contribution, take the 
form: 



(y/2(p - 20 - 2 PN ) + m d (V2(2f' - p) + 2p N ) ) tr x 



( V2(2£ + p) - 2 PN ) + m d (2p N - V2(2£' + p)) ) a 

Thus using the corresponding <5-functions the r\ and r\ integrations can be done trivially. 

• The exchange term in the one pion contribution. 
By this we mean that one of the interacting particles participates in the pion (see fig. [5]) 
. Proceeding as above have: 

^/3/37r = J^j3 S (P2+P3 -P2 ~ P 3 J 6 (pi - p^j 5 (p 4 + p 4 ) CT 4 . (p 4 + P 4 J (63) 

Going into the Jacobi variables we find: 

1 , /2 Pjv - >/2 (Vfy + W +3(€-C' + p)Y 

X/ r> r> r> \x I \ / 

(27T) 



(64) 



We find it convenient to use the above S functions to obtain: 

Z' = l(V2 q + 2V3 V + 6Z + 3 P -V2 PN ), ^ = 29 " 3 ^ +2PJV 

Thus the £ and ry can be done trivially. Furthermore A-terms, appearing in the case of 
the pseudoscalar contribution, take the form: 

<J2 (m d (q + V2 {VSv + 3(£ + p)) - 4pjv) + m u (a + V&V - 3\/2$ + 2p N )) 

A 2 = ^ ^ i '— '- ^ ^ (65) 

Qm d m u 

. 0-3 (m u (q - 2 v / 6?7 + 2p N ) - 3 (q + y/2p) m d ) 

A 3 = — i i '- i '- '- (66) 

bm d m u 



TABLE I: The spin flavor matrix elements of the various spin operators encountered in this work. They 
are normalized to the matrix element of the nucleon spin. 



( k indicates the spin ranks) 
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pseudoscalar 


5%/2 
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ia3 x 04 


axial 
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direct 


^2 
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((Ji. 0-2)0-4 


direct 


h 
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[(0-1 X o- 2 )fci2 = 2; 04] k = 1 


direct 


WlO 
9\/3 


(0-1.0-4)0-2 


direct 


_ 7V2 


(0-2.0-4)0-1 


direct 




0-4 


exchange 


J 


(0-2.0-3)0-4 


exchange 


* 


[(ct 2 x o- 3 )fc 2 3 = 2; 04] k = 1 


exchange 


SVlO 
9V3 


(0-2.0-4)0-3 


exchange 


13^2 
9 


(0-3.0-4)0-2 


exchange 


13V2 
3 



The Qvf3f3 decay amplitude at the nucleon level. 

Performing the orbital integrals we encountered in the previous section, we must evaluate the 
spin-flavor ME for the various operators encountered above, classified according to their spin rank. 
The obtained matrix elements, in units of the nucleon spin ME are included in table [J) . Using 
these results one can obtain the needed amplitude at the nucleon level. As expected from the above 
discussion we will consider three possibilities: 

1. The Of qq case 

In this case we can write the amplitude as 



1 1 ,/2m 

M = — m —, = VN.d ME(sf)J orb (67) 

where erjv is the nucleon spin and ME(s — f) is the spin- flavor matrix element ( see table [TJ) 
and J or b is the radial integral. One finds 



The coefficients Cj can be read off from Eqs 15111531 namely 



/ (3qm d + m u (2p w - 3q))\ p N ,„ n , 
^S(v) = 7. => (69) 

V bm,dm u J tun 



Cp= f (3qm M - m d (3q + 2pjy)) \ ^ jJ/v 



V 6mdm u J tun 



q a= ( (3q {m d + mj - 2m M pjy) \ ^ 3q-pjy 



V 6m d m u / m N 



The term pn of the amplitude will lead to a non local effective operator in coordinate space. 



2. The 3 P qq case 

Double beta decay proceeds via two quarks in a state with isospin one, which is color anti- 
symmetric. So the two quarks must be in a spin one state. So there is no contribution in 
V-A theories, since the vector and the axial vector contributions are identical. For the scalar 
and pseudoscalar cases the needed couplings depend on the particle model assumed. In the 
R-parity violating SUSY the coupling is , e.g. |(?7 T + |?7ps) found in [llj |. ln our discussion 
we will not include such a model dependent coupling. We will distinguish the two possibilities: 

a) The direct term. 

In this case we can write the amplitude as 



m =^wM Ai - mj - (72) 

In the case of the scalar contribution we find from table Q] that 

Ai.A 2 = --^-q.ffiv (73) 

In the case of the pseudoscalar contribution (see Appendix) and in the local approximation 
Pat = we find: 

Ai.A 2 = - \ — '- (<Ti.<7 2 )g.<7 4 (74) 

o \\/ / ix z + bbjym^mu/ 

We expect this to be a good approximation. In any event it makes the operator tractable. 
The corresponding orbital integral is: 

J u - 2332 Mg) r -bla 2 /6 r -blpl((2^)/(3+2x 2 ))/6 _ 

8jV3 - b % p %U2 x *)/( 3 +2x*))/6 (75) 

10\/2> 2 (0)tt 3 / 2 V ; 

We not with satisfaction that any uncertainties in the pion w.f. have dropped out, at least if 
the non local term in the exponential are ignored. 

b) The exchange term. 

The amplitude takes the form: 



M = 



1 1 V2r 
(2tt) 3 / 2 3V3V272 



A2.A3J, 



3 J orb 



(76) 



Again there is no contribution in V-A theories, since the vector and the axial vector contribu- 
tions are identical. In the case of the scalar contribution we find from table [J that 



. A V2 
A 2 .A 3 = — q.crjv 

y 

In the case of the pseudoscalar contribution for the constituent quark masses we get: 



A2.A3 



, 320^ (7x 2 + 1) (56a; 2 + 3) 2 416^2 (588x 4 - 77a; 2 + 57) 



147 (28a; 2 + 3) m 



+ 



/V 



63 (28a; 2 + 3) 2 b 2 N m 2 N 



a N .q 



(77) 



(78) 



where x = jf-. Note the presence of the q 2 in the first term. This will lead to an operator with 



a different radial dependence, i.e. (x) 
orbital integral for the exchange term is: 

3 3 2 7 %/2 



V 2 i^(a;) (see Eq. (TIT)]) ). The corresponding 



Jo 



9r 



(3 + 28a; 2 )V3 + 28a; 2 (<^(0)) 2 



(°) „ ((q-p«) 2 /6)((4x 2 )/(3+28x 2 )) 



(79) 



or 

Jorb = gr 648^/3 (2, 2 + 3) 3/2 e _ fcM(q _ PN)V6)((4 , 2)/(3+28 , 2)) (g0) 



5tt 3 / 2 (28a; 2 + 3) 3/2 ^(0)m p ^K 
In this instance the obtained results depend on the pion w.f. at the origin (via x). 

RESULTS 

Our main results are the coefficients a 27r and a\ n , which multiply the standard nuclear matrix 
elements. We will not elaborate further on the new non local terms (at the nucleon level). 

The coupling coefficients 



Before presenting our results we should mention that in the elementary particle treatment 11 1 
one can write 

2 j 2 ( m T 



a2 * = W A 9rK Wv) m 

Obtained under the factorization approximation: 

< 7T+| Jpjp\ir- >= - < 7T+| J P \Q >< 0\j P \iT- >, <Q\Jp\tt- >=mlK (82) 
The parameter h v is given by 

K = «\/20.668 — — — (83) 
m d + m u 

Returning back to our approach we note that the non relativistic reduction is applicable in the con- 
stituent quark mass framework, m u — rrid « mjv/3. In this case the pseudoscalar term contribution 
becomes: 

a 27r = -0.0005 (for x=1.0) , a 27r = -0.05 (for x=0.4) 

We should compare this with the value obtained in V-A theory see Eq. (f2"4"|). using f% NN = 0.08 
and 6jv = 1-0 fm: 

a 27r = 0.013 (x = 1) and a 27r = 0.11 (x = 1.0) 

i.e. it is quite a bit smaller. It is also much smaller than the value 0.20 obtained in the elementary 
particle treatment [llj using current quark masses. This disagreement cannot be healed by the 
fact that in the present case we encounter a very strong dependence of the results on the pion size 
parameter, see Fig. ??, unless we use very unrealistic values of the pion size parameter. One expects, 
of course, an enhancement of the pseudoscalar contribution, if one uses the current quark masses, 
since they are assumed to be very small. Indeed this way for typical values x = 1, &at = 1 fm, 
md = 5 MeV and m u = 10 MeV we obtain a 27T = — 1.3 and a 27r (17i 7r ) = 0.08, which are very large. 
We should mention, however, that the validity of the non relativistic reduction at the quark level 
may be questionable in this case. 

The coupling coefficients ai n 

Before proceeding further we will briefly present how the coefficient ai n was obtained in the 
context of the elementary particle treatment 



pros 

Q 



The needed parameters were obtained using the factorization approximation one writes in the case 
of 1 — 7r mode 

5 

< p\jpJp\nTr> — - < p\Jp\n >< 0|Jp|7r - >, < p\ J P \n >= F P w 4.41 (85) 
o 



The matrix element < 0| Jp\n > was given above (see Eq. (|82|l). Thus these authors 11] find: 

a l7r = -4.4 x 1(T 2 (86) 

Returning to our approach these coefficients are obtained in the following procedure: First we 
write 

- 2 M = c l7rffr ^ (87) 
nip 2m p 

Then, ignoring the momentum dependence in the exponential, we get: 

1. Double beta decay only. 

From Eqs([69 ]) -(|71 |) we see that the only local contribution comes from the axial current. 

ci-K = -=== Ji7r( x ) (current masses) (88) 

9g r \/b 3 N m d m N m u 

cijr = — / _^_ - — / lir (x) (constituent masses) (89) 

with 



Sgry/b^mff 



ft{x)=x* /2 (90) 
Using ma = 5 MeV, m u = 10 MeV and g r = 13.4 we get 

citt = 1.6 fi n (x), (current masses ) (91) 
On the other hand for the constituent masses we find: 

Ci,r = 3.4 x 10~ 2 f^(x), (constituent masses ), (92) 

The corresponding coefficient that must multiply the nuclear matrix element is 

f 2 

airr = citt 2 (93) 
J A 

aiTr = 0.085 fx v {x) (current masses ),ai 7r = 1.8 x 10~ 3 fi n {x) (constituent masses ) (94) 

2. The direct term 

There is no contribution of the direct diagram if the non local terms are ignored. 

3. the exchange term 

• In the case of the current quark masses we get the standard term: 

c ln = 1.0 x 10 3 /iT(^) (95) 



In addition we have an operator which results from the term in the amplitude, which was 
cubic in q. Thus we factor out the q 2 m^ and absorb it in the effective transition operator. 
In the remaining coefficient we merely replace q 2 by m 2 .. Thus 

c l7r = 91.5g^ r (x) (96) 

Proceeding as above we get respectively: 

a l7r = 51 fST(x) or a l7r = 51 g<% r (x) (97) 

The coefficient f^ r (x) is associated with the standard operator fli„(x„), while g^ r {x) 
must be linked with a new type of operator ^li„(x„) with modified radial dependence 
, i.e. F^ k \x) ==>• —V 2 F^ k \x) (sec Eq. (fTU|) ). Both coefficients are so normalized that 
/f" r (l)=<7^ r (l) = 1. In any case the use of current quark masses leads to very large 
values. 



• The constituent quark masses. 
In this case we get: 



c l7r = 1.37 AT(x) or c l7r = 1.72g%?(x) (98) 



a l7r - 0.071 fi° n (x) or a l7r = 0.090 5 f°"(x) (99) 

Again the coefficient /f°™(x) is associated with the standard operator, while gi° n (x) must 
be linked with the operator ^(i,), with f[° n {l) = g1° n (l) = 1 - 

The functions /£, / c " r l7r(x), firgf r (x) f^,f1° n {x) and ^°™(x)are shown in Fig. For i = 1 
for the standard local lit operator considering all contributions mentioned above with con- 
stituent quark masses we find a\-K = 7.3 x 10~ 2 , which is in size almost a factor of 2 larger 
than that obtained in elementary particle treatment [U (see Eq.®) . Note, however, that 
our results depend on the pion size parameter. 




FIG. 7: The functions which provide the dependence of 1 — 7r amplitude on the pion size parameter through 
the variable x = b^/bN are exhibited. On the left we show the relevant coefficients using the current quark 
masses. The continuous curve is associated with the coefficient ( see Eq. (|94[) h the long dash is associated 
with with the exchange q-independent coefficient (/i" r ) and the the short dash with that of gi^ r (see Eq. 
(|97|) ). On the right we show the same quantities obtained with constituent quark masses. 



DISCUSSION 



In the present paper we have considered the effective Qv(3(3 decay operator associated with the 
exchange of heavy particles mediated by pions in flight between nucleons. A harmonic oscillator 



non relativistic quark model in momentum space was employed for the pion and the nucleon. This 
allowed one to separate out the relative from the center of mass motion. The ratio of the pion to 
the nucleon harmonic oscillator parameter, x — b^/bjy was treated as a parameter. When needed, 
the constituent quark mass equal to 1 /3 of the nucleon mass employed. 

The obtained results were compared to the elementary particle treatment, with current quark masses, 
previously employed. 

In the case of the two pion mode we find a new term with different momentum dependence, which 
is not present in the elementary particle treatment. This gives rise to a new operator, which has the 
same structure as the one previously associated with the one pion mechanism. 

In connection with one pion mechanism we found that there exist three diagrams, which cannot 
be distinguished in the elementary particle treatment, namely: 

1. Diagrams in which the qq is crated out of the vacuum via the strong interaction. 

In this case we employed the 3 Po model. The strength of this interaction was fitted to the 

pion nucleon coupling g r . 

We distinguished two possibilities: 

• The two interacting quarks participate only in the structure of the nucleon. 

• One of the interacting quarks participates in the structure of the pion. 

2. The the qq is crated by the weak interaction itself. 

Depending on the mechanism we encountered new non local terms, i.e. terms which depend on the 
nucleon momentum. These will lead to new types of effective nuclear operators, which have not 
been examined up to now. 

The results obtained in the present calculation depend among other things on the ratio of the pion 
to nucleon size parameters. Using reasonable values for this ratio we obtain values of which 
are in good agreement with those obtained in the elementary particle treatment. Regarding the 
couplings a2-K, however, we find that they are slightly smaller than those obtained in the elementary 
particle treatment in the case of the V-A theory. They are, however, quite a bit smaller than those 
obtained in the case of the pseudoscalar term, when the constituent quark masses are used. We 
can, of course, obtain much larger values for the pseudoscalar term, if the current quark masses 
are used. Admittedly, however, it may not be very consistent to do so in our approach, since it is 
essentially a non relativistic treatment. We thus suspect that the small current quark masses are 
behind the large values found in the elementary particle treatment. 

In summary, taking into account the fact that a number of approximations are behind both 
approaches, we may say that there exists a reasonable agreement between them, which gives a 
degree of confidence in both. A more complete comparison can, of course, be made only after the 
inclusion in the calculation of the nuclear matrix elements of the new operators found in the present 
approach, namely: i) the local operator J7i 7r (x 7r ) resulting from terms cubic in q and ii) the non local 
operators, which depend on the nucleon momentum. 
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APPENDIX 

In this appendix we will present the relevant formulas in the case the qq pair is produced in via 
the strong interaction. In the case of the 1-pion contribution we get: 



1. The direct term 

In this case we can write the amplitude as 



M = TT^E^^^ A ^ J ^ (100) 
(2^)3/2 373 



In the case of the scalar contribution we find from table Q] that 

/2 

Ai.A 2 - --^-q.ajv (101) 
In the case of the pseudoscalar contribution we find: 

A A {x 2 + 2) (m d - m u ) PN-&1 (x 2 + 2) (m d - m u ) PN-&2 ( 2p N \ ^ 
(2x 2 + 3) m d m u (2x 2 + 3) m d m u \ 2x 2 + 3 / 4 

1 (x 2 + 2) (m d - m u )p N <ri (m d -m u ) 2^/2 

-02. 04 



3 (2a; 2 + 3) 

nid m u \[2\j2x 2 + 3bj\!m d m u y/2x 2 + 36/v 

1 (x 2 + 2) (m d - m u )p N <r 2 (m d -m u ) 2^2 

3 (2a; 2 + 3) mdTn,u \/~2\j2x 2 + 3b^m d m u \/2x 2 + 36jv 



CT1.CT4 



1/ 2piv \ / (m d -m u ) \ 2 
+ 3 [ q + 2x^3) * [ V2V2X 2 + 3b N m d m J (1 ° 2) 

In the limit = m u the above expression vanishes. In the local approximation pjy = we 
find 

a a 1 ( (m d -m u ) \ 3 

Ai.A 2 = - v i (en. (72)9.04 (103) 

o \v4a; / + 6bNm d m u J 

We expect this to be a good approximation. In any event it makes the operator tractable. The 
corresponding orbital integral is given by Eq. (175p . 

2. The exchange term 

The amplitude takes the form: 



m =^w 3 "tB am - (104 » 

Again there is no contribution in V-A theories, since the vector and the axial vector contribu- 
tions are identical. In the case of the scalar contribution we find from U that 

V2 

A 2 A 3 = — q.a N (105) 

y 

In the case of the pseudoscalar contribution we find: 

(7a; 2 + 1) (21 (14a; 2 + l) m d - (70a; 2 + 9) m„) ((224 x 2 - 9)m d + 3 (224a; 2 + 19) m„) 



A2.AC, = -- 



441 (28a; 2 + 3) J m 2 d m 2 u 
g.cr 2 g.CT 3 g.CT4 

((1400a; 2 + 117) m d - 3 (168a; 2 + 23) m u ) ((448a; 2 + 45) m u - 2lm d ) 

3528 (28a; 2 + 3) 3 m 2 m 2 u 
4 (llTO d + bm u ) (21 (14a; 2 + l) m d - (70a; 2 + 9) m u ) 



PN-O-lPN-O-iPN-O^ 



63 (28a; 2 + 3) 2 b 2 N m 2 rn 2 u 



a 2 .(Ji q.cr 3 



2 (1 \m d + 5m tt ) (21ra d - (448a 2 + 45 ) m u ) 
63 (28a; 2 + 3) 2 b 2 N m 2 m 2 u 



2q.<j 2 (7m d + m u ) ((224a; 2 - 9) m d + 3 (224a; 2 + 19) m u ) 



63 (28a; 2 + 3) 2 b 2 N m 2 m 2 u 
2 (7m d + m u ) ((1400a; 2 + 117) m d - 3 (168a; 2 + 23) m„) 

o O3.O4 Vn.<J 2 

63 (28a; 2 + 3) 2 6^m 2 m 2 
8 (7a; 2 + l) (llm 2 + (4a; 2 + 7) m u m d + 2 (6a; 2 + l) to 2 ) a 2 a z 

t <72-<73 <7-<74 

3 (28x 2 + 3) 2 b 2 N m 2 d ml 
2 Pn .<t 4 (llm 2 + (4a; 2 + 7) m u m d + 2 (6a; 2 + l) m 2 



-0-2.0-3 PN-Vi (106) 



3 (28a; 2 + 3) 2 b 2 N m 2 m 2 u 
In the limit of ignoring the non local terms we get: 

(7a; 2 + 1) (21 (14a; 2 + l) m d - (70a; 2 + 9) m u ) ((224 a; 2 - 9) m d + 3 (224a; 2 + 19) m„) 



A2.A3 



441 (28a; 2 + 3) 3 m 2 m 2 

q.a 2 q.a 3 q. 04 

4 (llm d + 5m u ) (21 (14a; 2 + l) m d - (70a; 2 + 9) m u ) 



63 (28a; 2 + 3) 2 b 2 N m 2 d m 2 u 



-0-2-0-4 9.0-3 



0-3.0-4 q.a 2 



2q.a 2 {7m d + m u ) ((224a; 2 - 9) m d + 3 (224a; 2 + 19) m v 

63 (28a; 2 + 3) 2 b 2 N m 2 ml 

8 (7a; 2 + l) (llm 2 + (4a; 2 + 7) m u m d + 2 (6a; 2 + l) to 2 ) a 2 cr 3 

t 0-2.0-3 0.0-4 (107) 

3 (28x 2 + 3) 2 b 2 N m 2 m 2 u 

where x = In the special case m u — m d — we get for the local terms: 

64 (7a; 2 + l) (56a; 2 + 3) 2 256 (56a; 2 + 3) 
A2.A3 = ■ - — { — — q.a 2 q.a 3 q.o± + _ — ± ^T7^ Q- a 3 



441 (28a; 2 + 3) J m 2 N 63 (28a; 2 + 3)" b N , f . x 



256 (56a; 2 + 3) 32 (4a; 2 + 5) (7a; 2 + l) 



0-3. 04 g.02 —k -cr 2 .a 3 q.oA (108) 



63 (28a; 2 + 3) 2 b 2 N m 2 N 3 (28a; 2 + 3) 2 b 2 N m N 

while the non local terms become: 

16 (56a; 2 + 3) 2 256 (56a; 2 + 3 



2 

A 2 A 3 = ^ — - — 7 Pn-o-iPn-OzPn-Oa \ N 2 ,/ 9 a 2^A PN-&3 

% 63 (28a; 2 + 3) 2 b 2 N m 2 N 



256 (56a; 2 + 3) 8 (4a; 2 + 5) , , 

0-3.0-4 PN-&2 H ; -o— n ^-O- 2 .0- 3 pjv.O-4 (109) 



63 (28a; 2 + 3) 2 b 2 N m 2 N 3 (28a; 2 + 3) 2 b 2 N ,„ N 

The first term of the local equation can be cast in the more suitable form by noting that: 

cr 2 . 00-3. 00-4-0 = y ^0-2.030-4.0 ^2 [(o- 2 x o- 3 )fci2 = 2 x o- 4 ] fc_1 .q^j (110) 

Using the spin matrix elements of table U we finally get using the current quark masses 
A2.A3 = 5V2q 2 

(7x 2 + 1) (21 (14a; 2 + l) m d - (70a; 2 + 9) m u ) ((224a; 2 - 9) m d + 3 (224a; 2 + 19) m u ) 

1323 (28a; 2 + 3) 3 m 2 m 2 

+26V2 

(-7 (224a; 2 - 75) to 2 + 2 (1176a; 4 - 938a; 2 + 93) m u m d + (7056a; 4 + 2212a; 2 + 201) to 2 ) 

567 (28a; 2 + 3) 2 b 2 N m 2 d m 2 u 
a N .q (111) 



where is the nucleoli spin, while for the constituent quark masses we get 



Ao.Aa 



, 320^ (7a; 2 + l) (56a; 2 + 3) 2 416^ (588a; 4 - 77a; 2 + 57) 



147 (28a; 2 + 3) d m 



N 



63 (28a; 2 + 3f b 2 N m 2 N 



a N .q 



(112) 



Note the presence of the q 2 in the first term. This will lead to an operator with a different 
radial dependence, i.e. F^ k) (x) =*> -X7 2 Fj; k} (x) (see Eq. (fTI)))), 
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